A Novel Symmetric Four Dimensional Polytope Found Using Optimization Strategies 
Inspired by Thomson's Problem of Charges on a Sphere 
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Inspired by, and using methods of optimization derived from classical three dimensional elec- 
trostatics, we note a novel beautiful symmetric four dimensional polytope we have found with 80 
vertices. We also describe how the method used to find this symmetric polytope, and related 
methods can potentially be used to find good examples for the kissing and packing problems in D 
dimensions. 



Using an optimization method (described below) in- 
spired by ones ourselves and others have used (see 
UlQ lJEIIj] and refs. therein) for a problem in three di- 
mensional electrostatics-Thomson's 6] problem of find- 
ing the minimum energy of N unit point charges on the 
surface of a unit conducting sphere-we have found a novel 
beautiful symmetric configuration with 80 vertices (Fig. 
1). The polytope has 64 vertices with 12 nearest neigh- 
bors of 0.7624, 0.6707, 0.7654 and 0.6661 distances, 16 
vertices with 10 neighbors of 0.7654 and 0.6661 distances. 
We have not seen this polytope previously 0, H| 

We found the configuration looking at the slightly ar- 
tificial, in four dimensions, but potentially useful (see 
below), problem of finding the minimum energy con- 
figuration of N charges (points) on the surface of the 
hypersphere (S3) x 2 
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mensions with the energy function 2j=yV 
is the Euclidean distance between two points 1 a nd 2 
\/ (xi - x 2 ) 2 + (yi - V2) 2 + (zi - z 2 ) 2 + (wi - w 2 ) 2 . As 
most surely like in three dimensions [2j, in four or more 
dimensions the number of good local mimima for this 
problem grows exponentially with N, and thus we can- 
not be certain that for N = 80 or other N that we have 
found the global minima. Nevertheless, even good local 
minima may be interesting or important configurations. 
In this initial work we have found our best local mini- 
mum for a given N by starting from 100 random initial 
starting configurations and then used a standard conju- 
gate gradient optimization. We have looked at N = 2 to 
200. The other As for which we have thus far found nice 
symmetric configurations are for N = 5,8, 24 and 120 for 
which we found as our minimum energy configuration the 
simplex (4D equivalent of the tetrahedron), the 16 cell 
(4D equivalent of the octahedron), the 24 cell and the 
600 cell (4D equivalent of the icosahedron) , respectively 
(Fig. 2), four of the six completely regular Platonic solids 
in four dimensions. We did not find the other two regu- 
lar polytopes, i.e. N = 16 the tesseract (or hypercube, 
4D equivalent of the cube) and N = 600 the 120 cell 
(4D equivalent of the dodecahedron). Their geometries 





FIG. 1: A novel symmetric four dimensional polytope with 80 
vertices. Four views showing each a different type of nearest 
neighbor bond (different distance). Red dots correspond to a 
vertex with 10 nearest neighbors. These figures show a simple 
parallel projection horn 4D to 2D, i.e. z and w coordinates 
are discarded, then each point is plotted accordingly to its x 
and y coordinates. All of the views have been rotated to show 
a symmetric pattern in 2D. 



are not energy minima, similar to what happens with 
the cube and the dodecahedron in 3D Thomson's prob- 
lem (Ref. Q and references therein). Using a method 
related to ours other higher dimensional polytopes have 
been found j^]. 

The D dimensional kissing problem asks the number of 
non-intersecting D dimensional unit hyperspheres that 
can touch a central unit hypersphere of dimension D. 
The maximal number of such spheres is known as the 
kissing number for that dimension, ATd. The kissing 
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Simplex, pentatope or pentachoron Cross Polytope, 16 cell or hexadecachoron 




24 cell or icositetrachoron 600 Cell or hexacosichoron 



FIG. 2: Regular polytopes found with Thomson's problem 
minimization process. These figures are simple parallel pro- 
jections from 4D to 2D. 

number is only proven in a few dimensions: K\=2, K2 — 
6, K 3 = 12, Kg = 240, K 24 = 196569[lil. (#4 = 24 
seems also to have been proven pj|). See |l0j and refs. 
therein for the best known kissing number in other di- 
mensions. The approach of minimizing l/ r u i n 
D dimensions potentially may be useful in finding new 
larger kissing numbers. Indeed, if for a given D with 
largest known kissing number Ku it can be found that 
a minimum energy configuration with Kjj + 1 (or more) 
charges/points can be found with each pair of points at 
least D-dimensional Euclidean distance one from each 
other, then a new kissing number for that dimension is 
found. There is freedom in choosing the energy function 
to be minimized, e.g. taking E = J^jjj V r 5 where n is 
an integer > 1. Also we note that for a configuration of 
N points found in D dimensions with greatest separation 
s between pairs of points, one thus immediately has an 
analogue of a kissing problem of putting hyperspheres of 
radius at least s around a central unit hypersphere. 

Now, a hypersphere has positive curvature and thus 
is not completely the same as flat space. To use min- 



imization of electrostatic problems to study the D- 
dimensional packing problem-to find the highest density 
of D-dimensional hyperspheres that can be packed into 
an infinite D-dimcnsional (flat) space an analogue of the 
above procedure can be used by working on the D dimen- 
sional surface of a D + 1 dimensional torus: Indeed, for 
example, the smallest distance between points on the two 
dimensional surface of the standard, simple three dimen- 
sional torus (doughnut/bagel) determines the density to 
which two dimensional spheres/balls-i.e., disks-can be 
packed. 

The clearest parameterization of tori we have found for 
this situation is to describe the D dimensional surface of 
a (D+ 1 dimensional) torus by D coordinates {x\, Xn) 
where each Xi lies in [0, 1). (So each xi is an angle scaled 
to lie in [0, 1).) For any of the Xi coordinates define 
\\xi - yi\\ = \xi - yi\ if \x 4 - yi\ < 1/2 or 1 - \x 4 - yt\ if 
\%i — Ui\ > 1/2- Then the appropriate squared intrinsic 
distance between two points on the surface of the torus 
is ||(xi,x 2 , ...,x n ) - (2/1,2/2, ■■■,2m)H 2 = T,i ~ 2/i 1 1 2 - 
One then seeks configurations of points/unit charges that 
minimize the energy function E = Ylijtj V \\ x — U\\ 
where n is an integer > 1. After the minimization is 
complete define r = 1/2 • rain ||x — j/||, where the mini- 
mum is taken over all pairs of points. The volume of the 
torus is simply 1=1. The volume of the D dimen- 
sional hyperspheres which can be packed on the surface 
of the torus (and also the packing density/packing frac- 
tion since the torus volume equals 1) is N -Q,d -r D , where 
fl D = tt d / 2 -T {{D + 2)/2) and N is the number of points 
on the surface of the torus. From the packing of points 
on the D dimensional surface of the D + 1 dimensional 
torus one then immediately gets a packing for infinite D 
dimensional flat space by placing a hypersphere of radius 
r at every point in R D which differs from the coordinates 
of a point on the torus by integer amounts in each coor- 
dinate. 
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